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ABSTRACT 

We construct a large data set of global structural parameters for 1300 field and cluster spiral galaxies 
and explore the joint distribution of luminosity L, optical rotation velocity V, and disk size R at I- 
' and 2MASS X-bands. We focus our study on scaling relations and residual correlations of benefit 

, to galaxy structure and formation studies, rather than deriving precise distance estimators. The I- 

• and iiT-band velocity-luminosity (VL) relations have log-slopes of 0.29 and 0.27, respectively with 

o'lniVL) ~ 0.13, and show a small dependence on color and morphological type in the sense that 
redder, early-type disk galaxies rotate faster than bluer, later type disk galaxies for most luminosities. 
The VL relation at /- and X-bands is independent of surface brightness, size and light concentration. 
I The modeling of the size-luminosity (RL) relation is somewhat sensitive to the adopted fitting method, 

due in part to strong dependences on surface brightness and color and significant uncertainties in the 
measurement of disk scale lengths. The log-slope of the /- and iC-band RL relations is a strong 
I function of morphology and varies from 0.25 to 0.5 with a mean of 0.32 for all Hubble types. At 

. most luminosities, early-type disk galaxies have shorter scale lengths than later-type ones - this latter 

■ trend may be reversed at high luminosities. The average dispersion ain{RL) decreases from 0.33 at 

/-band to 0.29 at K, likely due to the 2MASS selection bias against lower surface brightness galaxies. 
Measurement uncertainties are CTinv — 0.09, crinL — 0.14 and somewhat larger and harder to estimate 
for R. We show that the color dependence of the VL relation is consistent with expectations from 
simple stellar population synthesis models; the RL relation shows a weak but opposite trend expected 
from those models, suggesting that at a given stellar mass, smaller galaxies are redder. The VL and 
RL residuals are largely uncorrelated with each other with a correlation coefficient r = —0.16 and 
Alog^jL / Alogi?|L = — 0.07 ± 0.01; the RV — RL residuals show a weak positive correlation with 
r = 0.53. These correlations suggest that scatter in luminosity is not a significant source of the scatter 
in the VL and RL relations. The observed scaling relation slopes, zero-points, and scatters can be 
understood in the context of a model of disk galaxies embedded in dark matter halos that invokes 
low mean spin parameters and dark halo expansion, as we describe in our companion paper (Button 
et al. 2007). We discuss in two appendices various pitfalls of standard analytical derivations of galaxy 
scaling relations, including the TuUy-Fisher relation with different slopes. Our galaxy data base is 
available at www.astro.queensu.ca/~courteau/data/VRL2007.dat. 

Subject headings: galaxies: dynamics — galaxies: formation — galaxies: kinematics — galaxies: spirals 
— galaxies: structure — dark matter 
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1. INTRODUCTION 

Understanding the origin and nature of galaxy scal- 
ing relations is a fundamental quest of any successful 
theory of galaxy formation. The success of a particu- 
lar theory will be judged by its ability to predict the 
slope, scatter, and zero-point of any robust galaxy scal- 
ing relation at any particular wavelength. Some observed 
scaling relations in spiral galaxies, based on their size, lu- 
minosity, and rotation speed, can be reproduced individ- 
ually to fairly good accuracy by invoking galaxy forma- 
tion models that include virial equilibrium after dissipa- 
tional collapse of spherical cold dark matter (CDM) ha- 
los and angular momentum conservation (e.g. Mo, Mao, 
& White 1998, hereafter MMW98; van den Bosch 1998, 
2000, hereafter collectively as vdBOG; Navarro & Stein- 
metz 2000, hereafter NSOO; Firmani & Avila-Reese 2000, 
hereafter FAROO; Shen, Mo, Shu 2002). 

One of the most firmly established empirical scaling re- 
lations of disk galaxies is the Tully-Fisher relation (TFR; 
TuUy & Fisher 1977); a tight correlation between the to- 
tal luminosity and the rotation speed of a disk galaxy. 
We here refer to the TFR as the velocity-luminosity 
(VL) relation. To date, no single CDM-based model of 
galaxy formation can simultaneously reproduce the slope, 
zero-point, scatter, and color trends of the VL relation, 
match the shape and normalization of the local luminos- 
ity function, and explain the sizes, colors, and metallic- 
ity of disk galaxies (see, e.g., vdBOO; Bell et al. 2003a; 
Button et al. 2007, hereafter D07). In addition, simul- 
taneously accounting for the mass and angular momen- 
tum distribution of spiral galaxies in a gas dynamical 
context remains a major challenge for hierarchical for- 
mation models (Navarro & White 1994; van den Bosch 
et al. 2002b; Governato et al. 2007). A complete theory 
of galaxy scaling relations awaits a fuller understanding 
of structure forming mechanisms and evolutionary pro- 
cesses (e.g. star formation, merging, feedback, and cool- 
ing prescriptions) in galaxies. Likewise, the fine-tuning 
of these galaxy formation and evolutionary models de- 
mands a careful examination of empirical scaling rela- 
tions of galaxies. 

In order to set up a framework for the study of galaxy 
scaling relations, we examine the correlations of param- 
eters related by the virial theorem, V'^ oc M/R. We con- 
sider three fundamental observables for each disk galaxy: 
the total luminosity L, the stellar scale length i?d of the 
exponential disk, and the observed circular velocity V. 
The stellar mass, M*, can be estimated from the luminos- 
ity by assuming a stellar mass-to-light ratio, T* = M*/L. 
The size-luminosity (RL) relation of galaxy disks is also 
expressed as L oc So^^j, where Eo is the disk central 
surface brightness. 

Key to mapping fundamental dynamical trends in spi- 
ral galaxies, the measurement of VL and RL relations 
and detection of their correlated residuals require veloc- 
ity amplitudes measured at a suitably chosen radius rep- 
resentative of the flat part of resolved rotation curves, 
red/infrared magnitudes to minimize extinction and stel- 
lar population effects, accurate disk scale lengths, and, 
ideally, color terms from digital imaging with a broad 
baseline (e.g. B — K) to test for variations in the stel- 
lar population and extinction effects. One of the goals of 
this paper is to assemble such a data base for field and 



cluster galaxies (within the limits of available material). 

A study of scaling relations in irregular and spiral 
galaxies by Salpeter & Hoffman (1996) yielded the cor- 
relations Lb oc i?2-68 oc y^-^^ oc M^f^ a Ml^^, where 
Mhi, Mdyn, and R, are the Hi and dynamical masses, 
and a characteristic radius, respectively. The blue lumi- 
nosities, as used in that study, are notoriously sensitive 
to dust extinction and stellar population plus dynamical 
effects cannot be simply isolated. A new study of scaling 
relations in the (near-) infrared would provide more ro- 
bust dynamical constraints to galaxy formation models. 
While the VL relation has been examined at nearly all 
optical-IR wavelengths (e.g. Strauss & Willick 1995; Ver- 
heijen 2001, hereafter VOl; Masters et al. 2006; Pizagno 
et al. 2007, hereafter P07), comparatively few analyses 
of the RL relation of spiral galaxies have been reported 
so far (Salpeter & Hoffman 1996; de Jong & Lacey 2000, 
hereafter dJLOO; Shen et al. 2003, hereafter ShOS). This 
is partly because the accurate disk scale lengths needed 
to calibrate the RL relation, at any wavelength, have 
only recently become available for large databases (e.g. 
Courteau 1996; Dale et al. 1999; MacArthur et al. 2003, 
hereafter MCH03). Furthermore, most investigations of 
the VL and RL relations have used grossly incomplete 
databases owing to the nature of the sample (e.g. single 
cluster) or selection limits (magnitude, surface bright- 
ness, diameter, etc.) The availability of light profiles 
for very large galaxy databases (e.g. Sloan Digital Sky 
Survey, York et al. 2000 [hereafter SDSS] and the Two 
Micron All-Sky Survey, Skrutskie et al. 1997 [hereafter 
2MASS]) heralds a new era for the study of galaxy scal- 
ing relations (e.g. Shen et al. 2003; P07) with tractable 
selection biases. 

The empirical VL relation is expressed as 

LocV, (1) 

with the near-IR log-slope a ~ 3.0 ± 0.2 (Willick 
et al. 1997; GiovaneUi et al. 1997, hereafter G97; 
Courteau 1997; hereafter C97; Courteau et al. 2000; VOl; 
Masters et al. 2006). Reported values of the log-slope a 
range from 2.8 in the blue to 4.0 in the infrared (Willick 
et al. 1997; TuUy & Pierce 2000; VOl; P07) for both high 
and low surface brightness galaxies (Zwaan et al. 1995; 
vol). However, log-slopes steeper than a ^ 3.5 in the in- 
frared typically result from small samples and excessive 
pruning on the basis of idealized morphology or kinemat- 
ics, a narrow range of inclinations, redshift cutoffs, etc. 
(Bernstein et al. 1994; VOl; Kannappan, Fabricant, & 
Pranx 2002, hereafter KFF02). The slope, scatter, and 
zero-point of blue VL relations are predominantly dom- 
inated by stellar population and dust extinction effects 
(e.g. Aaronson & Mould 1983; G97; WiUick et al. 1997; 
TuUy & Pierce 2000; P07) and on the techniques used 
to recover the major observables and fitting for funda- 
mental relations (e.g. Strauss & WiUick 1995; C97; VOl; 
Bell & de Jong 2001; KFF02). Because we are mainly 
interested in masses, rather than luminosities, we do not 
concern ourselves with the VL and other scaling relations 
measured at blue; wavelengths. 

It is shown in Appendix B that previous dynamical 
derivations of the observed VL relation which arrived at 
log slopes of 3 (—7.5 in magnitudes) based on the virial 
relation Mvir oc V^ (e.g., MMW98; NSOO) and 4 (-10 in 
magnitudes) based on disk dynamics used erroneous as- 
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sumptions. A complete physical interpretation of the VL 
relation is still missing (for a fuller discussion, see e.g., 
Gnedin et al. 2006; D07). In terms of observables, the 
most fundamental relation is between the total baryonic 
mass of a galaxy, inferred via its infrared luminosity and 
a stellar mass-to-light ratio and total gas mass (Hi -t- He), 
and its total mass inferred from the asymptotic circular 
velocity of the galaxy disk (McGaugh et al. 2000; Bell 
& de Jong 2001; vol; McGaugh 2005; Geha et al. 2006; 
Gurovich 2006; De Rijcke et al. 2007). The "baryonic" 
VL relation is expressed as 

Mbar« 1/'^'=". (2) 

Since disk gas mass fractions typically increase with 
decreasing luminosity (e.g., McGaugh & de Blok 1997), 
one typically has Obar > cl. The slope abar is sensitive to 
the method used to determine stellar mass-to-light ratios. 
Using stellar populations models one finds Obar — 3.5 
(Bell & de Jong 2001; McGaugh 2005), while adopting 
stellar mass-to-light ratios derived from MOND results, 
by construction, in Obar = 4 (e.g. McGaugh 2005). 

While the log-slope of the VL relation can be repro- 
duced fairly well by naive derivations from CDM-based 
structure formation models (e.g. MMW98; NSOO; D07; 
see also Appendix B), the predicted scatter can be large 
compared to the inferred "cosmic" scatter of ^ 0.25 mag 
in red/infrared bands (Willick et al. 1996; VOl; P07) and 
interpretations about its dependence differ (see below). 
Besides the basic understanding of the slopes of galaxy 
scaling relations, the dependence of their scatter has also 
been addressed by many, especially for the VL relation 
(e.g. Aaronson & Mould 1983; Giraud 1986; Rhee 1996; 
Willick et al. 1997; KFF02; Courteau et al. 2003), and 
can be used to set realistic constraints on structure for- 
mation models (Courteau & Rix 1999, hereafter CR99; 
Heavens & Jimenez 1999; FAROO; NSOO; VOl; Buchal- 
ter, Jimenez, & Kamionkowski 2001; Shen, Mo, & Shu 
2002; Gnedin et al. 2006; D07). While various trends in 
the scatter of the blue VL relation have been discussed 
in the past, correlations of the near-infrared VL resid- 
uals with inclination, size, concentration, gas fraction, 
or far infrared luminosity are only few and inconclusive 
(Aaronson & Mould 1983; VOl; P07). The dependence 
of the scatter in the near-IR relation on color and surface 
brightness is, however, still a matter of contention that 
we discuss in ^12. 31 The scatter in the RL relation is also 
addressed in ^2.3\ 

The study of scaling relations in galaxies has bene- 
fited from the two-pronged application of the VL rela- 
tion for the purposes of: (i) Estimating relative distances 
to measure deviations from the mean Hubble flow (see, 
e.g., Strauss & Willick 1995 and the reviews in Courteau, 
Strauss, & Willick 1999); and (ii) testing galaxy forma- 
tion and evolution models (Dalcanton, Spergel & Sum- 
mers 1997; MMW98; vdBOO; D07). The philosophy of 
sample selection and calibration differs in both cases. For 
cosmic flow analyses, the calibration and science sam- 
ples must be pruned mainly on the basis of morphology 
in order to minimize systematic errors and ensure the 
smallest possible magnitude (distance) error. Since scat- 
ter in the VL relation depends strongly on the slope of 
the relation, it is found that the combination of steepness, 
magnitude errors, extinction correction, and sky stability 
favors red (i?&/) bands for smallest distance errors and 



cosmic flows applications (C97; Giovanelli et al. 1999; 
VOl). The accuracy of bulk flow solutions also depends 
on the size of the sample. In order to collect large enough 
samples, VL calibrations for flow studies rely mostly on 
Hi line widths or Ha rotation curves that can be collected 
relatively quickly on modest aperture telescopes. These 
rotation measures based on radio and optical spectra typ- 
ically sample the disk rotation out to ~3 — 4 disk scale 
lengths and are readily calibrated on the same system 
(C97; Dale et al. 1999). 

In contrast, the use of the VL relation as a test bed 
for galaxy formation models requires the widest range 
of morphological types, to sample all structural proper- 
ties, and that extinction and stellar population effects be 
tractable to isolate genuine dynamical correlations. The 
nearly dust-insensitive infrared bands are thus ideal for 
such applications (VOl). Rotation velocities are prefer- 
ably extracted from fully resolved Hi rotation curves ob- 
tained using aperture synthesis maps that sample the 
disk rotation out to roughly 4—5 disk scale lengths. 

Ideally, the study of scaling relations should rely on ho- 
mogeneous complete volume-limited samples assembled 
with the very purpose of testing for broad structural and 
dynamical differences; at the moment, we must contend 
with the more finely pruned heterogeneous samples of 
early and mostly late-type spirals that have been col- 
lected during the last decade mostly for flow studies. It 
is currently not possible to correct any of the existing 
catalogs for incompleteness or assemble a volume-limited 
sample extending to low mass that is complete in any 
meaningful sense (e.g. McGaugh 2005). The latest het- 
erogeneous data bases, which may include near-infrared 
luminosities and colors, for large samples of spiral galax- 
ies still enable us to characterize the dependence of scal- 
ing relations under the assumption that the different se- 
lection biases even out when sufficiently many different 
samples are collected together, and to examine various 
constraints of structure formation models in ways hith- 
erto unsettled. 

Courteau et al. (2003) showed that barred and un- 
barred galaxies have similar dynamical properties (be- 
yond the co-rotation radius) and that they share the 
same VL relation. As a natural extension of this study 
and CR99, we use in §2.21 and §2.31 extensive all-sky 
distance-redshift catalogs presented in fj^l to character- 
ize the mean VRL disk galaxy relations. In an attempt 
to extend our analysis to the (nearly) dust-free domain, 
we consider in §2.4l the infrared photometry from 2MASS 
for a large sample of spiral galaxies and discuss any rel- 
evant caveats. Although the 2MASS misses low lumi- 
nosity and low surface brightness galaxies it still serves 
some purposes for our analysis of scaling relations. In 
^2.5\ we bolster the notion of surface brightness and size 
independence of the VL relation, and we address in ^ 
the scatter of the VL and RL relations and their weak 
dependence on disk color. The main results are summa- 
rized in 21 Four appendices give A) figures for the mean 
VRL relations for each data sample, B) a derivation on 
the origin of disk galaxy scaling relations, C) an alter- 
native derivation of the TuUy-Fisher relation, and finally 
D) color transformations. 

The collected data and basic observations in this pa- 
per form the basis of the analytical modeling of spi- 
ral galaxy scaling relations in our companion paper 
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(D07. Our basic data base is available publicly at 
www.astro.queensu.ca/~courteau/data/VRL2007.dat. 

2. AVAILABLE DATA AND BASIC ANALYSIS 

We consider four major samples for which accurate 
near-infrared galaxy observables and rotational velocities 
are available. These include (a) the large /-band survey 
of galaxy distances for bright field spirals in the Southern 
sky by Mathewson, Ford, & Buchhorn (1992; hereafter 
"MAT" ) ; the all-sky /-band VL surveys of cluster and 
field spirals by (b) Dale et al. (1999; hereafter "SCII"), 

(c) Courteau et al. (2000; hereafter "Shellflow"); and 

(d) , the multi-band BRIK VL survey of Ursa Major 
cluster galaxies by TuUy et al. (1996) and VOl (hereafter 
"UMa"). The first three VL surveys were originally de- 
signed to map the convergence of the velocity field within 
~60ft.~^ Mpc, with the SCII and Shellflow studies pay- 
ing special attention to VL calibration errors for data 
from different observatories. By design, these surveys fa- 
vor bright late-type galaxies with inclinations on the sky 
greater than 45°; most have i ~ 60°. 

Properties for each sample are given in Table 1. These 
include in Col. (2), the number of galaxies from the orig- 
inal database that have the full complement of relevant 
observables; Col. (3), the nature of the sampled galax- 
ies (cluster or field; the predominantly "field" MAT and 
Shellflow surveys include a small fraction of cluster galax- 
ies); Col. (4), the original photometric band coverage. 
i?-band magnitudes for all the MAT galaxies were ex- 
tracted from the RC3 (de Vaucouleurs et al. 1991). We 
also have SDSS g — i colors for 39 SCII galaxies (see Ap- 
pendix C) and JHK 2MASS Kron magnitudes and col- 
ors for the brightest 360 SCII galaxies; Col. (5), the mag- 
nitude or diameter limits of the original catalog; Col. (6), 
the redshift limits of the survey in kms~^. For UMa, 
we adopted the cluster distance of 20.7 Mpc {HST Key 
Project, Sakai et al. 2000); and in Col. (7), the rotation 
measure, extracted from spatially resolved Hi or Ha ro- 
tation curves. 

While the MAT and SCII samples used both Hi 
linewidths from integrated spectra and Ha rotation ve- 
locities from resolved rotation curves, we use only the 
sub-samples of galaxies with available Ha, which is 
largest for both samples. For the purposes of galaxy 
structure studies Ha rotation velocities are prefered over 
Hi linewidths as the rotation velocity is measured at 
a known point of the rotation curve. Furthermore, Hi 
linewidths suffer from uncertain turbluence corrections. 

Rotation velocities for the Shellflow sample were mea- 
sured from resolved Ha rotation curves. For the UMa 
database, we use the sample of 38 galaxies with Vmax 
measured from resolved Hi rotation curves (VOl). 

Disk scale length measurements were computed in 
somewhat different ways for each sample. For the MAT 
galaxies, we have used the two-dimensional bulge-disk 
(B/D) decompositions of MAT /-band images by Byun 
(1992). These decompositions assumed a de Vaucouleurs 
bulge profile which is not ideal for late-type galaxies 
(MCH03) and therefore the extrapolated disk central sur- 
face brightnesses may be biased low by a few tenths of 
a magnitude. Those scale lengths are however deemed 
adequate for the purpose of our study since they match 
the mean RL relation (see ^2.2^ . SCII disk scale lengths 
were obtained by fitting a straight line to the exponen- 



tial part of the /-band surface brightness profile from 
~21 to ~25 /-mag arcsec"^ (the so-called "marking 
the disk" technique). Shellflow scale lengths were ex- 
tracted from one-dimensional B/D decompositions of 
azimuthally-averaged /-band surface brightness profiles 
(Courteau et al. 2000; MCH03). These fits account for a 
Sersic bulge and an exponential disk. Disk scale lengths 
for the UMa sample were measured from a "marking the 
disk" technique but the fit baseline is unspecified and 
erratic fits are reported in TuUy et al. (1996). The lat- 
ter could be due to inclusion of bulge light in the disk 
fit (scale lengths biased low) for the brighter galaxies 
and sky underestimate (scale lengths biased high) for the 
fainter galaxies. However, much like MAT, these pecu- 
liarities do not affect the UMa RL relations as compared 
against the other samples (see §2.21 and Appendix A) . 

The conversion of apparent observables to physical pa- 
rameters uses the distance scale adopted by the original 
authors, except for UMa galaxies for which we adopted 
D = 20.7 Mpc (Sakai et al. 2000). 

Deviations from a smooth Hubble flow may cause dis- 
tance errors, especially for the nearby galaxies in UMa. 
For the more distant galaxies in MAT, Shellflow, and 
SCII, typical thermal velocities of ~200 kms^^ would 
cause relative distance error of ~3%. We find in H2.1\ 
that distance errors which affect mainly the luminosi- 
ties and sizes are small compared to the other sources of 
observational uncertainty. For more details about each 
sample, the reader is referred to the original papers. 

2.1. Corrections 

In order to homogenize the four data samples, we have 
applied the inclination corrections for the line widths, 
luminosities, and scale lengths in a uniform way, as de- 
scribed below. 

The rotation velocity estimates are corrected for incli- 
nation and redshift broadening using 

V = Vobs/[{smi)il + z)]. (3) 
The inclination angle is derived via 

where the semi-minor and semi-major axes are deter- 
mined from CCD isophotal fitting, and go is the intrinsic 
flattening ratio of edge-on spiral galaxies, here taken to 
be go = 0.2 (e.g. Haynes & GiovaneUi 1984; Courteau 
1996; TuUy et al. 1998; Sakai et al. 2000). 

The apparent (/-band) magnitudes, to/, are corrected 
for both internal extinction, Ai^t, and external (i.e.. 
Galactic) extinction, ^ext- A small k-correction, A]^, is 
also applied such that 

mi = TO/,obs - ^int - ^cxt - ^k- (5) 

Internal extinctions are computed using the line width 
{W = 2V) dependent relation from TuUy et al. (1998): 

^int = liiW) log(a/6) = [0.92 + 1.63(logiy - 2.5)] log(a/5) 

(6) 

where a/b is the major-to-minor axis ratio and ^i(W) is 
a line width dependent extinction parameter. The ex- 
ternal (Galactic) extinction is computed using the dust 



Scaling Relations of Spiral Galaxies 



5 



extinction maps of Schlegel et al. (1998), while the k- 
corrections are computed using the line width dependent 
formalism of Willick et al. (1997). 
The absolute magnitudes, M, are computed using 

Af = m-51ogi?L-25; D^^ = [1 + z) , (7) 

with VcMB the systemic velocity of the galaxy in the 
reference frame at rest with the cosmic microwave back- 
ground (Kogut et al. 1993). The /- and if -band lumi- 
nosities are computed using for the Sun Mj^q = 4.19 and 
Mk^q = 3.33. 

Disk scale lengths from /-band images are corrected 
for inclination using 

i?i = i?i,obs/[l + 0.41og(a/6)], (8) 

(Giovanelli et al. 1994), and converted into kpc using the 
distance Dl- The correction above (Eq. ([8])) is somewhat 
uncertain but the choice of observed or deprojected scale 
lengths does not affect our final results. 

Finally, central surface brightnesses are corrected for 
inclination. Galactic extinction, and cosmological dim- 
ming (per unit frequency interval) using 

Mo,/ = Mo,/,obs + 0.51og(a/5)-Aoxt-2.51og(l-fz)^ (9) 

The coefficient 0.5 in front of the extinction term log(a/5) 
is empirically determined by demanding that the resid- 
uals of the relation between the /-band central surface 
brightness and rotation velocity has no inclination de- 
pendence^. We have assumed full transparency in the 
/f-band. Following Giovanelli et al. (1997) we assume 
an uncertainty of 15% in 7/, a/6, A^xt, and ^k, and 
propagate the errors. Including distance uncertainties, 
the average errors on the observables are (T\nv — 0.09, 
(Tini — 0.14, and ainR = 0.15 (adopting a fitting uncer- 
tainty of 15% on disk scale length). 

Distance errors in the VL relation distance affect 
mainly luminosity. If we include a distance uncertainty 
in terms of the recession velocity, the TF scatter is dom- 
inated by the uncertainty on V , so that linear fits with 
errors on V only or orthogonal fits with errors on both 
V and L are almost identical. Distance errors for the 
RL relation can be of issue since they affect both R and 
L. However, since R (x D and L oc D^, distance er- 
rors move RL data points along lines of constant surface 
brightness. Since the RL log-slope is ~0.3, the effect of 
distance errors is minimised. This is especially true given 
the larger intrinsic scatter and systematic uncertainties 
in determining scale lengths and correcting for internal 
extinction. 

Hubble types for all but the SCII galaxies were ob- 
tained from the heterogeneous NASA Extragalactic Data 
(NED) base. Most SCII galaxies are not classified in 
NED and their morphological types were determined in 
a homogeneous fashion via a combination of eye-ball ex- 
amination, B/D ratio and/or concentration index (Dale 
et al. 1999). Interacting and disturbed galaxies were re- 
jected in all samples. 

In addition to the VRL parameters and Hubble types, 
we have also considered for our study based on optical 
images the fully corrected V — I color, /-band central 
surface brightness, ^o,/, and concentration index. 

^ For a disk of zero thickness, one expects this factor to range 
from for an optically thick disk to 2.5 for an optically thin disk. 



For the Shellflow and UMa samples V — I and B ~ I 
colors were obtained from the original CCD photome- 
try. For the SCII sample V — I colors were derived from 
the SDSS (see Appendix C). For the MAT sample B-I 
colors were derived from RC3 B magnitudes. We con- 
verted B — I colors into V — I colors using the relations 
between color and stellar mass-to-light ratio in Bell & 
de Jong (2001). These yielded a total of 742 galaxies 
with available V — I. The central surface brightness /io./ 
are inward extrapolations of the disk fits described above 
and are available for all galaxies in our compilation. 

The galaxy light concentration index can be defined 
as: 

C72^rj5/r25, (10) 

where the radii enclose 75% and 25% of the total light ex- 
trapolated to infinity. For reference, a pure exponential 
disk has C72 =2.8, regardless of its total mass or scale 
length. Concentration indices for the MAT sample were 
not available. In total, 691 galaxies in our compilation 
have C72 estimates. 

Ahogether, the MAT, SCII, Shellflow, and UMa sam- 
ples combine for a total of 1303 separate entries. The 
range of physical parameters for our combined data set 
is shown in Fig. [T] This range is broad in every respect, 
yet not complete in any statistical sense. We have not en- 
deavored to correct each sample for incompleteness due 
to the complex inherent selection biases. Most parame- 
ter distributions can be closely approximated by a log- 
normal function. The relations between morphological 
type and luminosity, rotation velocity, size, color, sur- 
face brightness and concentration are shown in Fig. [2l 
The distributions for all the parameters, but /io,/, are 
broadly similar with a gentle decline from the highest 
mean values for the earliest-type systems (Sa) down to 
intermediate types (Sbc), settling onto a constant mean 
level for types Sc and later. This global trend is roughly 
reversed for /io,/. Our study of galaxy scaling relations 
will be extended to the /f-band based on 2MASS data 
in i j2.4l for comparison purposes. 

2.2. Mean Parameter Relations 

We consider the joint projected distribution in each 
of the planes defined by a pair of the three log virial 
variables, log/Zd, logV, and log Lj for the four combined 
samples in Fig. [31 the distributions for the individual 
SCII, MAT, Shellflow, and UMa samples are shown in 
Appendix A. Different symbols identify the full range of 
spiral Hubble types. 

The fits shown as solid black lines in each VRL figure 
are the orthogonal linear fits to the full combined data 
sets of 1303 galaxies. Table 2 gives the results of the or- 
thogonal flts using propagated errors on both variables 
(see e.g. Akritas & Bershady 1996) for the full sam- 
ple and for each Hubble type; the errors on the slopes 
and zero-points correspond to 1-a deviations assessed by 
bootstrap resampling. The regressions were performed 
over the full available range of luminosities, sizes, and 
velocities. The number of galaxies per Hubble type is 
shown in parentheses. The need for more measurements 
of the earliest and latest type spirals is obvious. 

We choose orthogonal flts over forward, inverse, or bi- 
sector fits for the following reasons. Forward and in- 
verse fits make the assumption that scatter exists only 
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in the dependent variable. Since it is not possible for 
both forward and inverse fits to be correct, bisector fits, 
which average the forward and inverse fits, can also not 
be correct (indeed, the bisector fit of a perfectly uncorre- 
lated distribution of two variables has an absolute slope 
of 1.) Furthermore, unlike orthogonal fits, forward and 
inverse fits cannot account for measurement uncertain- 
ties on both variables simultaneously. 

The combined /-band data give a VL relation, V oc L", 
with log-slope aj = 0.29 ±0. Of. The agreement between 
the log-slopes for the VL relations for different samples, 
as judged from Figs.fTSHTSlin Appendix A is fairly good. 
Our VL fit is also a close match to those published by the 
original authors (as Mj oc logT^), even though the min- 
imization techniques can be quite different. The Pear- 
son correlation coefficient of the VL relation is r = 0.92 
and the conditional distribution of F|i is roughly log- 
normal. Our VL relation, which uses optical rotational 
velocities, is a close match to the /-band VL relation of 
Masters et al. (2006) for their SF1++ catalogue of 807 
cluster galaxies when allowance is made for line width 
differences^ . 

The correlations for the RL and size-velocity (RV) 
relations are weaker, averaging The RL re- 

lation for the full sample, R (x L"", has a log-slope 
Pi = 0.32 ± O.Of . Since the fits of the VL and RL rela- 
tions are more robust than the RV relation, we adopt an 
RV relation requiring self-consistency between the three 
relations. That is, ii V (x L" and R oc L^^ and R (x V , 
then J = f3/a. 

Size-dependent distributions are less robust than the 
VL relation owing to such factors as, (a) lack of a uni- 
form and universal definition of scale length (MCffOS), 

(b) a morphological dependence that we discuss in 

(c) larger intrinsic scatter due to the natural disper- 
sion of the spin parameter A (D07), and (d) suscepti- 
bility to surface brightness selection effects. Whether 
scale lengths are measured from B/D decompositions or 
"marked" over a specified range of surface brightnesses 
(while omitting the bulge region or not) can yield scale 
differences greater than 20% (MCH03). Again, the over- 
lap between the four samples is too small to assess sys- 
tematic errors. It is therefore difficult to determine how 
much of the lower correlation coefficients in the RL and 
RV relations is genuine and/or due to inadequate scale 
lengths'^. Despite potential pathologies with the MAT 
scale lengths (see fj2]), the scaling relation slopes for this 
sample and SCII are comfortably close. In spite of in- 
trinsic differences between each sample (esp. selection 
functions), we make the assumption that our calibration 
in i 32.f [ brought all the samples on the same system. 

It has been proposed that the scaling relations may be 
different for the fainter, LSB galaxies (e.g., Kauffmann 
et al. 2003; Shen et al. 2003), where gas fractions are 

^ The 7-band VL relation of Masters et al. (2006) is based on 
HI line widths and normalised to Sc galaxies. Application of their 
transformation from Ha to HI (see their Eq. 2) results in a steeper 
slope for the VL relation based on HI line widths. This explains 
in part the difference between the shallower log-slope of our VL 
relation (0.29) and theirs (0.32). 

^ The improved reliability of scale length measurements must 
come from an extensive program of deep near-infrared imaging to 
measure homogeneous and accurate disk scale lengths and half- 
light radii for a very large collection of spiral galaxies based on a 
single reduction method and high signal-to-noisc data. 



likely to be different. Such a departure is not apparent 
in the samples that most closely probe the LSB regime 
(SCII and MAT; see Figs. [Hl&Iiniin Appendix A). Our 
observations are also corroborated by the study of low 
mass dwarf galaxies of Geha et al. (2006) who finds no 
break in the TF relations of very faint and brighter spiral 
galaxies. 

In summary, we find that the mean /-band scaling re- 
lations for spiral galaxies in Table 2 are: 

V^cxL°-29, /?dOcL°-32, ^^^^1.10 (11) 

The uncertainties in the log-slopes for the SCII sample, 
which we use as our baseline, are ±0.01 (VL), ±0.02 
(RL), and ±0.12 {RV) based on 1-ct measurement errors. 

A study of the volume-corrected RL distribution for 
140,000 late-type SDSS spirals (Shen et al. 2003; here- 
after Sh03) has yielded log-slopes from P ~ 0.23 (low 
L) to P = 0.53 (high L), with scatter ranging from 
c\nR — 0.45 (low L) to 0.30 (high L). These results are 
fully consistent with ours; see Table 2 where our /-band 
RL log-slopes agree perfectly with those of Sh03. This 
nearly perfect agreement must be contrasted in light of 
various caveats including the fact that the size measure- 
ment in Sh03 uses the Sersic half light radius, fit to a 
circularized SB profile. The median half light radius of 
their sample is 2", compared to a mean seeing FWHM 
of 1'.'5 which is cause for concern. The morphological 
sampling of the SDSS study also differs significantly from 
that of our collected samples. The size measurements are 
thus different from ours and possibly called into question 
due to seeing contamination and poorly matched profile 
shapes (see e.g. MCH03). 

The MAT sample (see Fig. was also analysed by 
dJLOO who attempted to correct for statistical incom- 
pleteness; they found the range p = 0.2 — 0.3. Our 
(non-corrected) result for the MAT sample alone is P — 
0.28 ±0.01 consistent with the upper envelope of dJLOO. 
The completeness correction affects mostly the low L 
distribution. At a given luminosity, one magnitude of 
surface brightness corresponds to 0.2 dex in R^, and a 
(complete) population of LSB galaxies would give a shal- 
lower log-slope in agreement with dJLOO. The effect is 
subtle, given the large intrinsic size errors, but we must 
keep that caveat in mind. 

2.3. The Morphological/Color Dependence of Galaxy 
Scaling Relations 

A simple understanding of the morphological depen- 
dence of the VL and RL relations is provided by Fig. SI 
We detect a slight and a strong morphological depen- 
dence of the VL and RL relations, respectively. This is 
confirmed by examination of Table 2. The morpholog- 
ical dependence of the VL relation was first discussed 
by Roberts (1978), and revisited by Aaronson & Mould 
(1983), Rubin et al. (1985), Giraud (1986), Pierce & 
TuUy (1988), G97, and KFF02 to name a few. It goes in 
the sense of early-type spirals rotating faster than later 
ones at most optical luminosities. The strong morpholog- 
ical dependence of the RL relation is such that earlier- 
type spirals are more compact than the later ones at 
most low to intermediate luminosities. This trend is re- 
versed for the brighter galaxies with L > lO^^L©. For 
early type (Sa) galaxies, the RL log-slope is close to the 
nominal limit of 0.5 for constant surface brightness (since 
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So oc L/R^). RL log-slopes are progressively shallower 
from early- to late-type systems. 

We can refer to Fig. [5] to assess whether a correlation 
with morphological type is driven by any other galaxy pa- 
rameters. On the basis of smallest scatter, the strongest 
correlations with morphological type would be with color 
(see also Roberts & Haynes 1994; Fig. 5) and velocity. 
Given that color is linked to the mean stellar M/L ratio 
and star formation history (SFH) of the galaxy (Bell & 
de Jong 2000; Bell et al. 2003a; Portinari et al. 2004), the 
latter might thus be a source of scatter in the VL relation 
(e.g. Heavens & Jimenez 1999; D07). One expects the 
VL scatter in bluer bands to be more sensitive to con- 
tributions of instantaneous star formation activity, while 
at (near-)infrared bands the scatter from stellar popula- 
tion effects mainly reflects the convolved star formation 
histories. 

The analysis of the Nearby Field Galaxy Survey 
(Jansen et al. 2000) by KFF02 identified B-i? color and 
Ha equivalent width as the main drivers of the scatter in 
the optical VL relation. Color and Ha equivalent width 
are both tributary of star formation histories though the 
former depends both on the integrated and instantaneous 
star formation while the latter is a function of the cur- 
rent star formation rate alone. We do detect a weak 
color dependence of the /-band VL relation (Figs. [5] & 
[5]). This result is further confirmed in Sj3]and Fig.[T31 and 
independently by P07. The trends are such that redder 
disks rotate faster and are more compact than bluer ones. 
However, while color may contribute to the VL scatter, 
we shall conclude in fJ3]that it is not a dominant source 
of scatter to either the VL or RL relations. 

Other than observational errors and scatter in lumi- 
nosity due to variations in the SFH, VL scatter is also 
expected from variations in the halo concentration pa- 
rameter and, to a lesser extent, scatter in the galaxy 
spin parameter and disk mass fraction. A discussion of 
the relative contributions to the VL scatter from SFH 
and halo parameters is presented in D07. 

According to VOl, all correlations of VL residuals 
with galaxy observables detected in blue bands, includ- 
ing color which affects the zero-point, vanish in the in- 
frared. This important assertion however rests on the 
study of a small data sample and deserves closer atten- 
tion. Contrary to findings by Aaronson & Mould (1983) 
and VOl, Rubin et al. (1985) found that the morphologi- 
cal type dependence of the VL relation was reduced, but 
not eliminated, at infrared H luminosities. In the section 
below, we make similar tests for any variations of the VL 
and RL relations at infrared wavelengths. 

2.4. The Infrared (2MASS) Velocity-Luminosity and 
Size-Luminosity Relations 

The advent of large-scale infrared surveys such as 
2MASS provides us with JHK luminosities, effective (i.e. 
half light) radii, and near-infrared colors to construct 
(nearly) dust-free scaling relations. However, with a typ- 
ical surface brightness limit of if ~ 20 mag arcsec"^, 
the 2MASS luminosity profiles are a full two magnitudes 
shallower than the typical J-band surface brightness pro- 
files for our galaxies. Still, the 2MASS data yield scaling 
relations that are as tight as the ones derived at /-band, 
with log-slopes representative of higher surface bright- 
ness systems. So while 2MASS data may be ill-suited 



for deep imaging investigations of galaxies (see e.g. Bell 
et al. 2003a) and deriving the correct slopes of the VL 
and RL relations (see, however, Kudrya et al. 2003), 
the information that they provide about the scatter of 
(2MASS) IR scaling relations may still be useful. Be- 
cause of the 2MASS magnitude limit of 13.5 /f-mag, only 
the 360 brightest SCII galaxies have measured 2MASS 
magnitudes. For simplicity, our analysis of /C-band scal- 
ing relations is restricted to the SCII sample. Note 
also that we make use of effective radii, supplied in the 
2MASS data pipeline, rather than disk scale lengths for 
this sample. For reference, a pure exponential disk has 
Re = 1.678/?d- 

We correct the IR data in a manner analogous to the 
optical data in ^12. II We use the iC-band correction from 
TuUy et al. (1998) scaled using the standard Galactic 
extinction curve (Schlegel et al. 1998) to get JHK ex- 
tinctions. These corrections remove any inclination de- 
pendences on the VL and color-L relations. We do not 
find any evidence of an inclination dependence to the 
values of i?e,A' and /ie,A' and thus no corrections to these 
quantities were applied. 

Figs. [9l-fTT] show the /C-band VRL relations for the 
360 brightest SCII galaxies. Based on this sample, the 
2MASS infrared scaling relations are V oc L^^'^, Re.K oc 
L°jf \ and Re,K oc V^-^^. Thus the /-band and 2MASS 
/C-band for the bright SCII sample have comparable VL 
relations, but the log-slopes for the RL and RV rela- 
tions differ substantially. We must keep in mind that 
the 2MASS RL and RV log-slopes are biased against the 
lower surface brightness SCII galaxies hence the steeper 
log-slopes for the RL and RV compared to the /-band. 
If we restrict the /-band sample of SCII galaxies to the 
bright 2MASS sub-sample, we find nearly identical /- 
band log-slopes for the VL (0.28) and R^L (0.36) rela- 
tions. This confirms that the different slopes for the /- 
and /sT-band RL relations result purely from selection, 
rather than a systematic difference between the /-band 
disk scale lengths and /sT-band effective radii, or putative 
pathologies with the 2MASS sizes. 

Because the VL relation is surface brightness indepen- 
dent, which we verify for 2MASS data in Fig. [Til we can 
still use these data to assess the parameter dependence 
in the VL and, to a lesser extent, the (incomplete) RL 
relation^. 

The /C-band VL and RL relations show similar de- 
pendences on morphological type as the /-band relations 
(Fig.[4|). As we see from Figs. l9l-fT2lthe slope, zero-point 
and scatter of the if-band VL relation is fully indepen- 
dent of surface brightness, size, and concentration pa- 
rameters. Any correlation with color, J ^ K this time, is 
difficult to assess given the paucity of late type systems 
in this sample; the V — I color dependence that was more 
prominent (though barely) in Fig. [8] seems to have van- 
ished at NIR wavelengths in Fig. [121 The RL relation 
shows the usual dependence on surface brightness and no 
dependence on J — /C color. 

* Although 2MASS surface brightness profiles only reach ~ 20 
mag arcsec"^, the 2MASS J — K color terms are still useful for the 
separation of redder and bluer disks since IR galaxy bulge and disk 
color gradients are already substantial after one disk scale length 
(MacArthur et al. 2004). 
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2.5. Surface Brightness Independence of the VL 
Relation 

Previous reports of the correlation of VL residuals 
with disk surface brightness may have generated confu- 
sion. Willick (1999) reported a correlation of /-band VL 
residuals with surface brightness and compactness for his 
"LPIOK" survey of distant cluster galaxies using a mo- 
ments fitting method to determine the exponential disk 
parameters even in the presence of irregularities in the 
galaxy light profiles. In an attempt to alleviate subjec- 
tive fitting boundaries, Willick's method used the entire 
surface brightness profile, including the bulge, thereby bi- 
asing the scale length and central surface brightness mea- 
surements. Applying his procedure, we can reproduce 
the putative surface brightness dependence of the VL 
relation while proper bulge-to-disk (B/D) fitting tech- 
niques (MCH03) find none. 

Prior to this study, the absence of surface brightness 
dependence of the VL relation had also been verified 
at optical bands by Sprayberry et al. (1995), Zwaan 
et al. (1995), and CR99, and in the infrared by VOL 
Our Figs. [6l \E\ fTTl and [T2l corroborate this evidence very 
clearly; namely that there is no dependence of the VL 
relation scatter on J-band central surface brightness /io,/ 
or if-band effective surface brightness. 

The surface brightness independence of the VL has 
been interpreted (e.g. Zwaan et al. 1995) as variation of 
the dynamical mass-to- light ratio, T, with surface bright- 
ness for a given total luminosity. If the surface brightness 
decreases, T increases in such a way that the VL rela- 
tion remains independent of surface brightness. A de- 
pendence of T on surface brightness is indeed expected 
in CDM-based galaxy structure models; the higher the 
disk surface brightness, the higher the contribution of 
baryons to the rotation velocity at, say, 2.2 disk scale 
lengths (Zavala et al. 2003; P05; D07). However, repro- 
ducing the surface brightness independence of the VL 
relation is a non-trivial task, especially for the highest 
surface brightness galaxies (CR99; D07). Indeed, for a 
given disk mass and dark halo profile, a small disk (higher 
surface brightness) would contribute significantly to V2.2, 
the total circular velocity at 2.2 disk scale lengths^. In 
the absence of halo contraction, the contribution of the 
halo to V2.2, for typical halo concentrations, slighty de- 
creases as disk size decreases. The net effect is typically 
that V2.2 increases as disk size decreases (CR99; DOT). 
When halo contraction is included a smaller disk results 
in a smaller halo and hence even higher total V2.2 (D07). 

Gnedin et al. (2006) circumvent this problem by as- 
suming that at a given stellar mass, smaller and thus 
higher surface density disks live in less massive halos. 
This solution is not particularly attractive as it has nei- 
ther observational nor theoretical motivation. On the 
other hand, the observed dependence of gas-to-stellar 
mass ratio on surface brightness is effective in reduc- 
ing the surface brightness dependence of the VL rela- 
tion (FAROO; van den Bosch 2000; D07). Furthermore, 
if galaxy formation involves clumpy cold flows, rather 
than smooth cooling flows, dynamical friction between 
the baryons and dark matter may reverse some (or all) 

^ The maximum velocity measurements collected in our data 
base correspond roughly to V'2.2- 



of the expected effects of halo contraction, thus reduc- 
ing the surface brightness dependence of the VL relation 
(D07). 

3. RESIDUAL CORRELATIONS OF SCALING RELATIONS 

In the spirit of CR99, we now examine below the corre- 
lations of the VL, RL, and RV residuals from the mean 
relations. We define the residuals for each object i as 
Ay{i) = y{i) — yfit(*)- The orthogonal fits for the global 
scaling relations on VRL were described in ^ 32.21 

Fig. [13] shows the residual correlations for combina- 
tions of AlogV, Alogi?, and AlogL for the four sam- 
ples based on /-band imaging (upper) and SCII sample 
based on 2MASS X-band imaging (lower). The corre- 
lation slopes and their associated errors are reported in 
Table 4 and are shown at the bottom of each panel in 
Fig. [131 The linear correlation coefficients, r, are also 
shown at the top left corner of each panel. The resid- 
ual correlations reported in Table 4 supersede those pre- 
sented in CR99 which differ slightly due to revised cor- 
rections to the raw observables. 

There is close agreement in the correlated residual so- 
lutions of all optical samples, with the general solution: 

A log V\L / A log i?|L = -0.07 ± 0.01 
Alogi?|y/AlogL|y= 0.71 ±0.02 
AlogF|i?/ AlogL|i?= 0.31±0.01. 

While we do not correct for possible sample incomplete- 
ness, the fact that all of our samples, with their very 
different selection functions, yield similar solutions sug- 
gests that any potential selection bias ought to be mini- 
mal or cancel out (see discussion in D07). The slopes for 
the AlogyjL — Alogi?|L residuals for all samples are 
negative (anti-corrclatcd) and low. This residual corre- 
lation is statistically different from A log V / A log R — 0, 
though only weakly. Note also that galaxies of all mor- 
phological types (and barredness, not shown here), scat- 
ter normally about the zero line (see Appendix A). 

FAROO have also computed the VL/RL residuals for 
the UMa sample (Fig. [T8|) and found no correlation for 
both for LSB and HSB galaxies, as do we for that sample 
(see Table 4 and the lower left panel of Fig. [T9)) . 

The weak surface brightness dependence of the VL re- 
lations has been used to infer the ratio of baryonic to 
dark mass in the inner parts of galaxies (CR99). D07 
however suggest that changes in the stellar mass-to-light 
ratio, variations of baryonic-to-stellar mass ratio, bulge 
formation prescriptions and assumptions about the re- 
sponse of the dark matter halo to the cooling baryons 
may all contribute to a weaker surface brightness depen- 
dence of the VL relation. 

With r ~ 0.5, the RV and LV residuals (middle col- 
umn of Fig. [13]) are also weakly correlated. The VR and 
LR residual distributions (third column of Fig. [T3[) are 
tightly correlated with r ~ 0.9 (a little less at K) and 
relatively small slope errors. This, as we explain below, 
is simply the VL relation recast in differential form. 

The correlations between residuals of the VRL rela- 
tions give clues to the contribution of scatter in V, L, 
and R. Consider a case where scatter exists only in V: 
then the Alogi?|y — AlogLjy residuals will be corre- 
lated with a slope equal to the log-slope of the RL rela- 
tion (i.e. ^^0.3). Similarly, for scatter only in L, the slope 
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of the AlogyjL — Alogi?|L residuals should equal the 
log-slope of the VR relation (i.e. ^1); for scatter only 
in R, the slope of the AlogT^ji? — AlogLji? residuals 
should equal the log-slope of the VL relation (i.e. ^0.3). 

Since there is intrinsic scatter in each of the VL, 
RL, and RV relations, there must be scatter in at 
least two V , L and R. The case of uncorrelated 
scatter in R and L, and no scatter in V, would re- 
sult in a positive A log V\L — A log R\L residual correla- 
tion, uncorrelated Alogi?|l^— A\ogL\V residuals, and 
a weak positive A\ogV\R — A\og L\R residual correla- 
tion. This situation is in clear conflict with all three 
observed residual correlations. The case with uncorre- 
lated scatter in V and R, and no scatter in L would 
result in uncorrelated Alog^jL — Alog_R|L residuals, 
A log A logLji? residuals with slope ^0.3, and pos- 

itively correlated A\ogR\V — AlogijV^ residuals. This 
latter is indeed in excellent agreement with the observa- 
tions. Thus in the simplest interpretation of the residual 
correlations, the scatter in the VL and RL relations is 
fully independent. However, as discussed in D07, repro- 
ducing independent VL and RL relations is non-trivial 
as the various theoretically expected sources of scatter in 
V\L and R\L often result in correlated V and R scatter. 

We conclude this section with a discussion of the color 
dependence of residual correlations. As discussed in ^2.3\ 
the VRL relations show a dependence on color and mor- 
phological type. Since both color and morphological type 
depend on luminosity and velocity, a most robust way to 
look for a color dependence of the VRL relations is to 
compare the VRL residuals with residuals of the color- 
luminosity and color-velocity relations. This is shown 
in Fig. [T31 The upper panels show the V — I color- 
luminosity and V ^ I color-rotation velocity relations. 
The middle panels show the VL residuals versus color 
residuals. The dashed lines show the predicted correla- 
tion slope for scatter in color at a given stellar mass, 
assuming log M^/Lj cx 1.26{V — /) from stellar pop- 
ulation models (Portinari ct al. 2004), and that scat- 
ter in color is uncorrelated with V and R. For the 
Alog V^|L/ — A{V — I)\Lj residuals, the observed corre- 
lation is weak but of the same sign as expected (from the 
stellar population models), in agreement with P07. This 
further confirms that scatter in color is a contributing, 
but not dominant, source of scatter in the VL relation. 
However, the A log Li\V — A{V — I)\V residuals are pre- 
fectly uncorrelated, also in agreement with P07. This is 
explained if scatter in the velocity-stellar mass relation 
depends on color, in the sense that redder galaxies of a 
given stellar mass rotate slower. However, this strikes 
against theoretical expectations that redder galaxies of a 
given total luminosity should rotate faster. 

We also find that the Alogi?|i/ - A{V - L)\Li resid- 
uals have only a weak correlation, but of opposite sign 
to that expected. This suggests that the scatter in the 
size-stellar mass relation depends on color, with redder 
galaxies of a given stellar mass being smaller. This is 
expected from theoretical galaxy formation models (Bell 
et al. 2003b). Finally the Alog R\V ~ A{V ~ I)\V residu- 
als show no correlation, but luminosity also does not play 
any role in this correlation, and therefore no dependence 
on colors is expected here. 

4. DISCUSSION 



The collection of some of the most extensive existing 
data bases of galaxy structural and dynamical parame- 
ters at /- and i^T-bands has enabled us to determine the 
following major observational results: 

• The /-band scaling relations of spiral galaxies, ob- 
tained by linear regressions of on i and R on L, 
are 

^^^o.29±o.oi^ i?d oc L?-32±o °i, i?dOcV^i i". 

• The 2MASS iiT-band scaling relations of bright 
SCH galaxies differ slightly (due to selection ef- 
fects) with 

T/ ^ r0-27±0.01 p ^ r0.35±0.02 n ^ t/1.29 

The slightly shallower log-slopes of the VL rela- 
tions from the /- to K-hand are expected based on 
the observed trend of color with luminosity. How- 
ever, the ii'-band RL and RV relations are biased 
to higher slopes due to the surface brightness bias 
of the 2MASS catalog for extended objects. A 
study of scaling relations based on deeper IR sur- 
vey data (e.g. UKIDSS; Hewett et al. 2006) should 
recover shallower slopes, comparable to those mea- 
sured here for /-band data, for the infrared RL and 
RV relations. 

• The log-slopes of the /- and if-band VL and RL re- 
lations show a dependence on morphological type, 
with steeper slopes for earlier type disk galaxies, 
especially so for the RL relation. 

• The /-band VL relation shows a dependence on 
color in the direction expected from simple stellar 
population synthesis models. However the RL rela- 
tion shows a trend (albeit weak) in the opposite di- 
rection, which is most simply interpreted as smaller 
galaxies at a given stellar mass being slightly red- 
der. 

• The VL relation shows no dependence on size or 
surface brightness at both / and K bands. 

• The AlogFji — A log/? I L residuals are weakly 
anti-correlated with correlation coefficient r ~ 
-0.2 and slope -0.07 ± 0.01. Unlike CR99, we 
showed in Dutton et al. (2007) that this result does 
not uniquely imply that most spiral disks have sub- 
maximal fraction of luminous to dark matter. 

• The Alog/?|F — AlogLjF residuals are weakly 
correlated with correlation coefficient r ~ 0.5 and 
slope ~ 0.7. Together with the weak Alog^li — 
Alog R\L residual correlation, this suggests that 
scatter in velocity and size dominate the scatter in 
the VL and RL relations respectively. 

• The AlogFj/?— Alog L\R residuals are strongly 
correlated with correlation coefficient r ~ 0.9 and 
slope ~ 0.3. This is the VL relation recast in dif- 
ferential form. 

The slopes zero points and scatter of the VL, RL 
and RV relations can be broadly understood in the con- 
text of disk galaxies embedded in dark matter haloes as 
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discussed in Appendix B and in more detail in Dutton 
et al. (2007). However, accurately reproducing the zero 
point and surface brightness independence of the VL re- 
lation, as well as the galaxy number density, may require 
a departure from the standard assumptions governing the 
formation of disk galaxies. 

In order to make further progress in constraining disk 
galaxy formation models from the VRL relations, a 
large {N > 1000) sample of galaxies with deep multi- 
wavelength imaging (near-IR is essential), high spatial 
and spectral resolution ID and 2D spectroscopy (with 
CO, Ha and Hi maps), controlled selection criteria, ac- 
curate and homogeneous measurements of structural pa- 
rameters, and stellar and baryonic masses, is needed. 
Ideally, chemical information should be available to con- 
strain stellar population effects. 
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APPENDIX 

/-BAND SCALING RELATIONS AND RESIDUAL CORRELATIONS FOR THE FOUR MAIN SAMPLES 

This appendix presents the separate distributions for the VRL scaling relations of the SCII, MAT, Shellflow, and 
UMa samples. These are shown respectively in Figs. [T5]fT51 The axis limits are the same in all figures and the fits 
shown are those for the full combined sample (see Fig. [3]). 

Fig- [m shows the residual correlations for combinations of Alog V^, Alogi?, and AlogL for the four samples based 
on /-band imaging. The colored point types have the same morphological dependence as in Fig. [3l All four samples 
give consistent residual correlations, namely a very weak AlogFji — Alogi?|L correlation (analagous to the surface 
brightness and size independence of the VL relation), a weak positive A log — A log L\V correlation, and a strong 
positive A\ogV\R — AlogLji? correlation (the differential VL relation). 
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ON THE ORIGIN OF DISK GALAXY SCALING RELATIONS 

The mean observed trends between V^bsj L and i?d at /-band are given by the following VL^ RL, and RV relations: 

Kbs(x£•'■'^ i?d«i°-^', i?d(xKibf. (Bl) 

In the spirit of the spherical-collapse model (Gunn & Gott 1972) one can define the virial radius, i?vir, of a collapsed, 
virialized gravitational body as the radius inside of which the average density is a factor Avir times the critical density 
of the Universe (the value of Avir depends on redshift and cosmology; see e.g., Eke, Navarro & Frenk 1998; Bryan & 
Norman 1998). Thus one can write the virial mass, Mvir, as 

Mvir = ^7ri?3i^Avi,.pcrit, (B2) 

where pcrit is the critical density of the Universe at a given z, and is given by 

where H{z) is the Hubble constant at redshift z. From the virial theorem, the circular velocity at the virial radius, 
Kir, is given by 

Kl = ^. (B4) 
Combining Eqs. lB2llB4l and setting H{z) = lOO/ikms^^ Mpc"^ we obtain 

M.„^fli„(.= (|5j)i (B5) 

/A ■ \'/' 

Kir = i?vir/lf^j , (B6) 



and 



-1/2 

M.ir = Vlh~'[^] G. (B7) 



Here G has units of (kms ^)^ kpcMg"'^, Afvir has units of A/q, i?vir has units of kpc, and Kir has units of kms ^. 

Thus, if A/vir/i, Rvir/Rd, and Vvir/Kbs are roughly constant for different galaxy types (i.e., not significantly corre- 
lated with either of the virial quantities), these virial relations are all that is required to explain the observed scaling 
relations of disk galaxies. However, there are many reasons to expect significant variation in Mvii/L, R^h/Rd, and 
Kir/Kbs- Below we realize these variations in terms of galaxy formation theory. 

The virial mass-to-light ratio can be written as 

A/vir/i = = — . (B8) 

/bar Cgf Csf "^d 

Here, T^, = M^,/L is defined as the stellar mass-to-light ratio, /bar is the baryonic mass fraction of the Universe 
(/bar = ^h/^M — 0.15 for a ACDM cosmology), Cgf is the galaxy formation efficiency that describes what fraction of 
the baryonic mass inside the halo ultimately ends up in the disk, and egf is the star formation efficiency that describes 
what fraction of the disk mass ends up in the form of stars. We combine /bar and Cgf into the disk mass fraction, md, 
and T* and egf into the disk mass to light ratio, Td- 

1 /3 

The virial relation Kir c>c A/^/. is thus related to the observed VL relation by the quantity 

^'obs 



CvL ^ ^ — . (B9) 

Kir V"^d/ 

In order to reproduce the observed VL relation Kbs oc L°'^^, Cvl must be proportional to Kbl"" ^^"^ has to reveal 
an amount of scatter that matches that of the observed VL relation. 

In the standard model for disk formation, set forth by Fall & Efstathiou (1980), Dalcanton et al. (1997) and MMW98, 
the relation between /?d and /?vir derives from the acquisition and conservation of (specific) angular momentum by 
the baryons, when cooling to form a centrifugally supported disk. The total angular momentum of a virialized system 
is conveniently expressed by the dimensionless spin parameter 

A = -^—hiL (BIO) 

(Bullock et al. 2001b). For an exponential disk embedded in a dark matter halo, the total angular momentum is given 

by 

Jd 27r / Y.[R)Vc{R)R^AR = A/d/?dKir/y (Bll) 
Jo 
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with and Vc{R) the surface density and circular velocity of the disk, and 

Jo •'vir 

(MMW98) a dimcnsionless number (e.g. /y = 2 if Vc{R) = V^ir). The specific angular momentum of the disk is related 
to that of the halo via 

If disk formation occured under conservation of specific angular momentum, as generally assumed, then Ad = A. 
Independent of this particular assumption one therefore obtains that 

^=V2^ (B14) 

-ftvir JV 

(MMW98). Thus, the virial relation i?vir oc V^ir is related to the observed RV relation by the quantity 

Chv ^ ^ f (B15) 



fv \ Vohs _ 

In order to reproduce the observed RV relation R^ oc V"^^ , Crv must be proportional to V^^~^, and has to reveal 
an amount of scatter that matches that of the observed RV relation. 

1 /3 

Similarly, the virial relation i?vir oc M^/. is related to the observed RL relation by the quantity 

Ad /Td^'/' 



Crl ^ ^ — . (B16) 
fv \mdj 

In order to reproduce the observed RL relation i?d oc L°'^'^, Crl must be proportional to L"rl-i/3^ g^^^^j ^las to reveal 
an amount of scatter that matches that in the RL relation observed. 

Thus, the challenge for galaxy formation theories is to understand how to meet the constraints on Cvl, Crvj and 
Crl derived here. That this is not a trivial matter becomes clear if one takes into consideration that (i) both md and 
Td depend extremely sensitively on the efficiencies of cooling, star formation and feedback, and are thus expected to 
be strongly correlated with Mvir and A (e.g., van den Bosch 2002a), and (ii) that both T^bs/Vvir and fv depend, in a 
convoluted way, on A, md, adiabatic contraction, and on initial halo concentration (e.g., D07), which itself depends on 
halo mass (e.g.. Bullock ct al. 2001a). 

In summary, while the virial relations of dark matter haloes provide a natural framework for understanding the 
observed scaling relations of disk galaxies, a direct linkage between the slopes of the observed and virial relations 
ignores the many physical processes that take place during galaxy formation. 
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AN ALTERNATIVE DERIVATION OF THE TULLY-FISHER RELATION 

Two simple predictions for the TuUy-Fisher {VL or TF) relation have been used in the past: one with a log-slope of 
3 (—7.5 in magnitudes), and the other with a log-slope of 4 (—10 in magnitudes). Here we compare both predictions 
and show that they are related to one another. 

From Appendix B, we have 

/A \ -1/2 / ^3 

^3 7 — 1 / ^vir \ / ^vi 



This is a TF relation with log-slope of 3, if Cvl is independent of V^bs (see also Dalcanton et al. 1997; MMW98; van 
den Bosch 1998; Syer, Mao, & Mo 1999). In contrast, several authors in the past have predicted that the TF relation 
should have a log-slope of 4 (e.g., Sargent et al. 1977; Aaronson, Huchra, & Mould 1979; Salucci, Frenk, & Persic 1993; 
Sprayberry et al. 1995; Zwaan et al. 1995). Their argument goes as follows: Assume that the disk surface brightness 
is described by an exponential function of the form 

/(r) =/oexp(-r/i?d), (C2) 
with i?d the disk scale length and Iq is the disk central surface brightness. The total luminosity of the disk is 

L = 2ttIoRI. (C3) 
Assuming that one measures the rotation velocity V^bs at a radius r ~ si?d, and that the gas and stars in the disk 
move on circular orbits. 



r 



6 



with M(r) the total mass within radius r. We can define the total mass-to-light ratio' 

T(0 . ^ . M,ir)+AUr)^ ^^^^ 

square (|C4p . and then use equations jCSj and (|C5|) to obtain: 

This is a TF relation with a log-slope of 4, as long as /q T^(r) is independent of Vohs and one measures the rotation 
velocity at a constant number of scale lengths^. The scatter is determined purely by the variation in Iq T^(r). This 
was emphasized by Zwaan et al. (1995), who, from the finding that both HSB and LSB spirals follow the same TF 
relation, concluded that LSB galaxies must have much larger values of T^(r) than their HSB counterparts. 

The TF relations in Eqs. (jCip and (|C6p thus predict a di fferent log -slo pe. T he cause of this apparent paradox lies 
in the fact that Iq and Vohs are not independent. From Eqs. IB5[ IB 141 and IC3| we derive the relation between surface 
brightness and observed rotation velocity 

Substituting this into Eq. IC6I yields a TF relation with a slope of 3 as long as 

^^""t^^/Fk;; ^^^^ 

is independent of Vohs- Rather than squaring Eq. IC41 a simpler expression results from substituting the relation 
between i?d and Vohs (see Appendix B) into Eq. IC41 thus resulting in 

\ , 1/ /A \ -1/2 
*\r\r I ^vir 



This is a TF relation with a slope of 3 as long as 

CvL^-f-^^ (CIO) 
T(r) Jv Vohs 

is independent of Vohs ■ 

Thus the two theoretical predictions of the TF relation, with log-slopes 3 and 4, are directly related. However, while 
Eqs. ICll and lC9l have an underlying TF slope of 3, variation of Cvl with Vohs (or L) will result in a steeper or shallower 
slope. For example, stellar mass-to- light ratios increase with Vohs (e.g. D07), and stellar mass fractions increase with 
Vohs (e.g. Kannappan 2004). This results in the steepest slope for the stellar mass TF relation, an intermediate slope 
for the baryonic mass TF relation, and the shallowest slope for the blue luminosity TF relation (e.g. Bell & de Jong 
2001). The slope of the TF relation will also depend on the velocity measurement, since the shape of galaxy rotation 
curves changes systematically with luminosity (e.g. C97; Catinella et al. 2006). 



® We use the subscripts "d" and "h" to refer to the disk and dark matter halo, respectively. 
Note that this derivation can be generalized to an arbitrary surface brightness profile by replacing the central surface brightness, Iq, 
with the effective surface brightness, 7(,ffi a^nd the disk scale length, R^, with the effective (i.e. half light) size, i?off- 
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SDSS COLORS FOR GALAXY CLASSIFICATION 

The need for dynamical and color information is crucial for any study of galaxy scaling relations. However, of all 
the available samples with homogeneous dynamical measurements (e.g. rotation curves) considered here, digital colors 
are only available for the Shellflow and UMa samples with VI and BRI imaging respectively. The MAT and SCII 
galaxy magnitudes were imaged only at /-band. In order to augment the available color base, we computed, in CR99, 
a B — I color term for MAT galaxies using B-band magnitudes for the RC3 (de Vaucouleurs et al. 1991). However, 
RC3 magnitudes are inherently uncertain with Am ~ 0.2 mag (Courteau 1996). 

We can extract SDSS g and i magnitudes (as proxy for V and / magnitudes) for the galaxies that overlap with 
our Northern samples. We do this for Shellflow and SCH galaxies in order to map the SDSS system onto the Lan- 
dolt/Cousins system upon which Shellflow and SCH magnitudes are based. We compare in Figs. [20lfc [211 instrumental 
Petrosian magnitudes from the SDSS fifth release against instrumental magnitudes for 31 Shellflow and 39 SCH galax- 
ies. The comparison of instrumental magnitudes avoids any differences due to extinction and k-corrections. Besides 
zero-point offsets, the g and i Petrosian magnitudes scale linearly with the V and / Cousins magnitudes in Shellflow. 
The same is true for SCH at near-infrared bands. Because the transformation slopes at / band are the same for 
Shellflow and SCH, we make the assumption that the transformation slopes at visual bands will also be similar. We 
also assume a monochromatic zero-point difference between the Shellflow and SCII magnitude systems of 0.13 mag. 

We compute transformations for the SDSS Petrosian g magnitudes into Cousins magnitudes using ^Pct — ^Land + 0.36 
for the SCII galaxies that overlap with the SDSS/DR5 database. The new instrumental "V" magnitudes from SDSS 
for 39 SCII galaxies were corrected the same way as the /-band magnitudes but using an extinction dependence 
Ay = 3.24, instead of Aj — 1.96, and the Poggianti (1997) formulation for the F-band k-correction. 
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TABLE 1 
Redshift-Distance Galaxy Surveys 



Sample 
(1) 


N 
(2) 


Gal. type 
(3) 


Phot, bands 

(4) 


Mag/Diam. limits 
(5) 


Redshift limits 
(6) 


Rot. measure 
(7) 


MAT 


545 


field 




Deso > 1-7 


< 7000 kms-i 


Ha 


sen 


468 


cluster 


I 


12 < mi < 17 


[5000-19,000] 


Ha 


ShelliJow 


252 


field 


V,I 


ruB < 14.5 


[4500,7000] 


Ha 


UMa 


38 


cluster 


B,R,I,K 


mz < 14.5 


20.7 ± 0.9 Mpc 


HI 



TABLE 2 

Orthogonal fits for the four combined /-band samples 



logy vs log L/ log Rd, I vs log L/ 



Type 


N 


slope 


zero-point 


"■in V|i 


slope 


zero-point 


'^In Ii \L 


ALL 


1303 


0.291 ±0.004 


-0.835 ±0.039 


0.132 


0.321 ±0.010 


-2.851 ±0.106 


0.325 


(V-I) 


742 


0.297 ±0.005 


-0.898 ± 0.056 


0.131 


0.306 ±0.013 


-2.681 ±0.136 


0.304 


Sa 


117 


0.303 ± 0.013 


-0.941 ±0.140 


0.126 


0.550 ±0.047 


-5.357 ±0.501 


0.375 


Sb 


570 


0.288 ± 0.008 


-0.805 ± 0.084 


0.132 


0.369 ± 0.020 


-3.373 ±0.213 


0.313 


So 


505 


0.272 ± 0.006 


-0.648 ± 0.060 


0.126 


0.328 ± 0.015 


-2.915 ±0.156 


0.320 


Sd 


111 


0.280 ±0.011 


-0.735 ±0.111 


0.125 


0.254 ± 0.021 


-2.116 ±0.213 


0.289 



TABLE 3 

Orthogonal fits for the SCH /C-band sample 









log V vs log Lk 






\ogRe,K vs log Lk 




Type 


N 


slope 


zero-point 


''"In V|L 


slope 


zero-point 


(^In R\L 


ALL 


360 


0.269 ±0.007 


-0.692 ± 0.081 


0.126 


0.346 ±0.018 


-3.13 ±0.19 


0.291 


Sa 


56 


0.274 ±0.021 


-0.721 ± 0.240 


0.129 


0.568 ±0.042 


-5.70 ±0.47 


0.266 


Sb 


131 


0.264 ±0.017 


-0.641 ± 0.185 


0.130 


0.428 ±0.035 


-4.06 ±0.39 


0.265 


Sc 


166 


0.255 ±0.012 


-0.549 ±0.124 


0.118 


0.388 ±0.025 


-3.54 ±0.27 


0.261 


Sd 


7 


0.249 ± 0.039 


-0.504 ± 0.404 


0.081 


0.314 ±0.168 


-2.68 ± 1.80 


0.358 



TABLE 4 
VLR Residual Correlations 







Alogy]L vs AlogR\L 


AlogK]V vs AlogL]V 


AlogVli? vs Alog L\R 


Sample 


N 


slope r 


slope r 


slope r 


ALL 


1303 


-0.07 ±0.01 (-0.16) 


0.71 ±0.02 (+0.53) 


0.31 + 0.01 (+0.93) 


MAT 


545 


-0.11 ±0.03 (-0.17) 


0.57 ± 0.04 (+0.59) 


0.32 + 0.01 (+0.90) 


sen 


468 


-0.09 ± 0.02 (-0.20) 


0.84 ± 0.08 (+0.53) 


0.32 ± 0.01 (+0.94) 


ShelWow 


252 


-0.08 ±0.03 (-0.18) 


0.78 ± 0.09 (+0.52) 


0.32 ± 0.01 (+0.93) 


UMa 


38 


-0.06 ±0.08 (-0.14) 


0.77 ± 0.17 (+0.48) 


0.31 ± 0.02 (+0.93) 


SCIIK 


360 


-0.03 ± 0.02 (-0.06) 


0.56 + 0.05 (+0.52) 


0.28 ± 0.01 (+0.88) 
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Courteau et al. 



Fig. 1. — Distribution of physical parameters for the four combined data sets. 



Fig. 2. — Correlation of galaxy observables (disk central surface brightness [^o = ~2.51og/o]i concentration C72, V — I color, luminosity 
Li, rotation speed V, and disk scale length Rd,i) with morphological type. The point types represent each samples: blue squares - Shellfiow, 
red stars - SCII, green circles - MAT, and black triangles for UMa. There is a weak dependence of all variables on Hubble type, with an 
apparent flattening of the trends for Sc galaxies and later. 



Fig. 3. — VRL scaling relations for all four samples color-coded by Hubble types. The linear orthogonal fits, shown by the solid lines, 
are reported in Table 2 for the combined samples; the 2-cf observed scatter are given by the dashed lines. 



Fig. 4. — VL and RL relations color-coded by Hubble types: (left) for the combined /-band sample, (right) for the SCH i^-band sample 
(the latter is discussed in i|2.4|l . The grey dashed lines show the 2-a scatter of the VL and RL relations for the full /— and i^— band 
samples. At 7-band, R = R^, the disk scale length; at /sT-band, we used R = Re, the eff'ective radius of the galaxy supplied by 2MASS. 
The VL relation shows, on average, that earlier-type spirals rotate faster than later-types at a given luminosity. The RL relation depends 
strongly on morphology with earlier types showing a steeper (^constant surface brightness) slope. 

Fig. 5. — Same as Fig.|4]but with point types representing different V — I color. The VL relation shows a weak dependence on color. 



Fig. 6. — Same as Fig. |4] but with point types representing different central surface brightness. The VL relation is fully independent of 
surface brightness, whereas the RL relation shows the trivial dependence on surface brightness. 



Fig. 7. — Same as Fig. |4] but with point types representing different concentration index. The VL relation is fully independent of 
concentration, whereas RL relation shows a weak dependence. 



Fig. 8. — Residuals from the VL (left) and RL (right) relations as a function of disk central surface brightness, disk scale length, 
rotation velocity, concentration index, V — I color, and optical morphological type for all the galaxies in our sample. Besides the expected 
dependence of the RL relation on /Jo, the scatters of the VL and RL relations may show weak correlations on color and concentration 
respectively. 



Fig. 9. — VRL scaling relations of the SCH K-hand sample separated by Hubble types. The linear orthogonal fits, shown by the solid 
lines, are reported in Table 3 for the combined samples; the 2-a observed scatter are given by the dashed lines. 



Fig. 10. — Same as Fig. [9]but with point types representing different J — K colors. The VL relation shows a weak dependence on color, 
while the RL and RV relations show none. 

FIGURE CAPTIONS 



N.B.: The full document, with B&W an d colou r figures, is available at: 

http:/ /www. astro. queensu.ca/~courteau/papers/VRL2007ApJ.pdf . 



Scaling Relations of Spiral Galaxies 
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Fig. 11. — Same as Fig. [9] but with point types representing different effective surface brightnesses. The VL relation is independent of 
surface brightness while the RL relation shows a trivial dependence on surface brightness. 



Fig. 12. — Residuals from the 2MASS VL (left) and RL (right) relations as a function of effective surface brightness fi^ft, effective size 
(/?eff)i rotation velocity, concentration index, infrared J — K color, and optical morphological type. The sample is grossly incomplete for 
types later than Sc. 



Fig. 13. — Correlations between the residuals of the VRL relations for the /-band sample (upper) and /C-band sample (lower). The best 
fitting slopes with l-cr uncertainty and correlation coefficients are given in the lower and upper left corner of each panel. There is only a 
weak correlation between V\L and R\L (left) compared to the prediction for pure disk only model (dashed line with slope —0.5). There is 
a significant correlation between the L\V and R\V residuals (middle), and even tighter correlation between V\R and L\R, which is simply 
the V L relation in differential form. 



Fig. 14. — Correlation residuals of the I-band relations with V — I color. The upper panels show the color-luminosity (left) and 
color-rotation velocity (right) relations. The middle panels show the residuals of the VL relation versus the color residuals at constant 
L (left) and constant V (right). The lower panels show the size residuals versus the color residuals at constant L (left) and constant V 
(right). The dashed lines give predictions based on scatter in color at fixed stellar mass and assuming, based on stellar population models, 
that color is related to stellar mass-to-light ratio. 



Fig. 15. — Scaling relations for SCII galaxies. Line widths are measured from Hq rotation curves and disk scale lengths are measured using 
the "marking the disk" technique (see text). Luminosities are computed from fully-corrected /-band magnitudes. Unlike Dale et al. the 
luminosities do not include any correction for morphological type dependence. The solid black lines and reported slopes correspond the 
orthogonal fits to the combined data set. The point types and colors correspond to the Hubble classification from the authors (Dale 
et al. 1999). 



Fig. 16. — Scaling relations for MAT galaxies. Line widths are measured from Ho rotation curves and disk scale lengths a re m easured 
from bulge-to-disk decompositions. Luminosities arc computed from fully-corrected /-band magnitudes. The fits are as in Fig. 1151 Hubble 
types are taken from the NASA Extragalactic Data Base. 



Fig. 17. — Scaling relations for Shellfiow galaxies. Line widths are measured from Ho rotation curves and disk scale lengths are measured 
fro m bu lge-to-disk decompositions. Luminosities are computed from fully-corrected /-band magnitudes. Fits and point types are as in 
Fig. [15] 



Fig. 18. — Scaling relations for UMa galaxies. Line widths are measured from HI synthesis rotation curves and disk scale lengths are 
measured using the "mar king the disk" technique (see text). Luminosities are computed from fully-corrected /-band magnitudes. Fits and 
point types are as in Fig. 1151 



Fig. 19. — Correlation residuals from the mean relation at fixed virial quantity for the Shellflow, SCII, MAT, and UMa samples. Symbols 
are as in Fig. [3] Galaxies of all shapes and types (barred and un- barred, early and later types) obey the same surface brightness independent 
scaling relations. The residual distribution for A log V / A log R of all three samples is flat for barred and un-barred galaxies of all Hubble 
types and luminosities. The middle figures show a weak correlation between L and R, and the figures on the right side recast the VL 
relation in its differential form. 



Fig. 20. — Comparison of raw SDSS (Petrosian) and Shellflow (Cousins) magnitudes. The two photometry data sets are derivable from 
one another, independent of a color term at least over the Shellfiow magnitude range. 



Fig. 21. — Comparison of raw SDSS (Petrosian) and SCII (Cousins) magnitudes. The linear transformation of magnitudes holds over 
nearly 3 magnitudes, independent of a color term. 



